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WATER  WAVES  PRODUCED  BY  EXPLOSIONS 
Herbert  C.  Kranzer 
and 
Joseph  B.  Keller 
Institute  of  Mathematical  Sciences 
New  York  University 

1,   Introduction. 

In  this  report  we  examine  the  production  of  water  waves 
by  the  following  mechanisms: 

1.  An  arbitrary  initial  distribution  of  impulse  applied  on 
the  surface . 

2.  An  initial  elevation  or  depression  of  the  surface  of 
arbitrary  shape . 

Later  we  will  consider  another  mechanism: 

3.  A  submerged  expanding  and  contracting  spherical  cavity 
(e.g.  an  underwater  explosion  bubble). 

These  mechanisms  are  intended  to  be  models  of  wave  pro- 
duction by  explosions.   The  first  model  applies  to  an  explosion 
above  the  surface,  the  second  to  a  surface  or  near-surface  explo- 
sion, and  the  third  to  a  deep  underwater  explosion.   In  the  sur- 
face case  the  explosion  will  produce  a  local  depression  of  the 
surface  and  waves  will  be  formed  as  the  water  flows  in  to  fill  the 
depression.   Since  this  model  neglects  the  initial  outward  flow  which 
accompanies  formation  of  the  depression,  the  model  with  an  initial 
elevation  might  be  more  realistic  because  it  yields  an  initial  out- 
ward flow.   However,  the  first  model  -  that  of  a  distribution  of 
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impulse  applied  on  the  surface  -  Is  perhaps  still  better,  even  for 
surface  explosions. 

In  our  analysis  we  assume  that  the  surface  elevation  h  , 
and  the  accompanying  water  motion,  are  symmetric  about  an  axis  .■ 
through  the  center  of  the  disturbance.   If  r  denotes  distance  from 
this  axis  and  t  denotes  time  then  we  may  write  h  =  Vi(r,t),   We  also 
assiirae  that  the  water  is  of  finite  xoniforra  depth,  which  we  denote 
by  h  ,  and  that  the  linear  theory  of  surface  waves  is  applicable. 
Both  of  these  assiimptions  will  be  analysed  in  a  later  report. 

In  the  next  two  sections  we  present  and  discuss  the  results 
of  our  analyses  for  the  cases  of  initial  impulse  and  initial  dis- 
placement respectively.   Typical  graphs  of  the  wave  height  as  a 
function  of  r  and  of  t  are  included.   Then  in  section  IV  these  re- 
sults are  applied  to  the  problem  of  wave  production  by  explosions. 
The  mathematical  analysis  is  contained  in  an  appendix. 

A  list  of  references  to  previous  work  on  this  subject  is 
appended.   Of  the  previous  work  the  earliest,  that  of  Poisson  [3], 
is  the  most  relevant.   However  he  considers  water  of  infinite  depth. 
Furthermore  he  obtains  detailed  results  only  for  paraboloidal  and 
ellipsoidal  depressions.   Burnside  [I4.]  also  considers  infinite  depth, 
but  only  twodimensional  motion.   He  obtains  a  series  for  the  wave 
motion  produced  by  an  initial  elevation  of  very  special  shape. 
Lamb  [$]    considers  only  the  case  of  infinite  depth.   He  presents  the 
point  source  solutions  for  initial  elevation  and  impulse  in  two  and 
three  dimensions.   Most  of  these  results  had  already  been  obtained 
by  Poisson  [3]»   The  contents  of  Lamb's  paper  are  reproduced  in  his 
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text  [1],  except  for  the  three  dimensional  results.   Rayleigh  [8] 
discusses  the  general  question  of  one  dimensional  propagation  of  a 
pulse  in  a  dispersive  medium.   Thompson  (Lord  Kelvin)  [12]  presents 
graphs  of  a  special  solution  describing  two  dimensional  waves  in 
water  of  infinite  depth.   (His  paper  [13]  is  a  reprint  of  [12].) 
In  [11]  he  discusses  the  method  of  stationay  phase  and  applies  it 
to  the  waves  produced  by  a  concentrated  elevation  in  water  of  in- 
finite depth.   Pidduck  [9]  obtains  a  power  series  for  the  solution 
resulting  from  a  concentrated  elevation,  and  discusses  instantaneous 
propagation.   In  [10]  he  presents  the  solution  for  waves  in  a  com- 
pressible fluid,  and  shows  that  instantaneous  propagation  results 
from  neglecting  compressibility.   Unoki  and  Nakano  [11;]  consider  waves 
due  to  an  initial  elevation  and  an  initial  impulse,  each  constant  over 
an  interval  of  the  surface.   They  treat  only  two  dimensions  and  in- 
finite depth.   Nevertheless,  their  results  are  in  fair  agreement  with 
observations  on  waves  caused  by  the  eruption  of  a  submarine  volcano. 
Terazawa  [?]  considers  the  motion  at  the  center  for  some  special  cases 
of  initial  impulse  and  initial  elevation  in  water  of  infinite  depth. 
He  treats  a  submerged  explosion  by  assuming  that  it  produces  an  initial 
distribution  of  impulse. 

On  the  basis  of  the  above  description  of  previous  work  on 
water  wave  production,  it  appears  that  the  production  of  three  di- 
mensional waves  in  water  of  finite  depth  has  not  been  treated.   There- 
fore we  have  treated  this  problem  for  arbitrary  initial  impulse  and 
initial  elevation,  using  the  method  of  stationary  phase. 
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2.    Initial  Distribution  of  Impulse  Applied  on  the  Surface. 

The  height  of  the  surface  following  an  Initial  impulsive 
force  distributed  over  it  is  given  by 

T  /R 

(2.1)     ^^(r,t)'~  -2 A  sin  2tt(|  -  ^)    for  r  »  R  . 

r/g  r 

In  (2.1),  I  denotes  the  value  of  the  initial  impulse  at  the  origin 
and  R  denotes  the  effective  radius  of  the  initial  impulse  (see  (2.9)].  The 
amplitude  A,  the  period  T  and  the  wavelength  \  of  the  waves  all  de- 
pend upon  the  ratio  r//gh  t  ,  and  A  also  depends  upon  the  initial 
shape  of  the  surface.   The  amplitude  A  is  zero  if  r/Zgh  t  >  1  , 
a  consequence  of  the  fact  that  the  maximum  group  velocity  of  any  wave 
is  v/gh 

Prom  (2.1)  we  see  that  the  surface  wave  height  decreases 
inversely  with  r.   The  wave  pattern  is  determined  by  the  slowly  vary- 
ing amplitude  factor  A  .   Each  portion  of  the  pattern  travels  outx^ard 
with  constant  speed  but  different  parts  travel  with  different  speeds. 
The  outer  parts  travel  faster  than  those  nearer  the  center.   Therefore 
the  pattern  is  lengthened  or  spread  out  as  it  travels  outward.   This 
lengthening  is  proportional  to  t  or  to  r.   Thus,  for  example,  the  dis- 
tance between  two  successive  maxima  or  minima  of  JAJ,  or  between  two 
successive  zeros  of  A,  increases  linearly  with  t  or  with  r  as  the 
waves  travel  out. 

The  amplitude  A  merely  yields  the  envelope  of  the  wave 
height.   The  factor  sin  2Tr  (^  -  ■^)    describes  the  actual  waves  -  an 
oscillation  of  height  which  is  rapid  compared  with  the  relatively 
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slow  rate  at  which  A  changes.   Thus  the  wave  height  is  given  by  a 
modulated  cosine  curve.   However,  the  period  and  wavelength  of  the 
waves  are  not  constant.   Both  the  wavelength  and  period  are  infinite 
at  the  outermost  part  of  the  wave  motion  (at  r  =  ^gh  t)  and  both 
decrease  as  r  decreases,  for  fixed  t,  or  as  t  increases,  for  fixed 
r.   This  variation  is  in  accordance  with  the  fact  that  waves  travel 
faster  the  longer  their  wavelength. 

The  functions  A,  T,  and  X   are  most  easily  expressed  in 
terms  of  the  auxiliary  variable  d  =   -^  ,   This  variable  is  itself 
a  function  of  rX/gh  t  and  is  defined  as  the  unique  non-negative  root 
of  the  equation 


(2.2)  ^(o)    =  |/i^     .      „„„   \,/2  ^ 


1/  ^\j/'^     slnha'  r^r  . 

2(cos   cf)-^'  /gh  t 


In  terms   of   <f  vie  have 


(2.3)  ^     =     ^  ' 


{2.k)  T      =  ^^ 


Vh 


tanho" 


1  /  0(tf)tanh  (f       ^  -,(fs 

o 


(2.5)  A  = rA-rro   /  '°i°^yn  °      tfKg),       r  < /IH  t  , 


0   ,  r   >  /gh  t 


.!  '■■% 
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The  function  !(-)  is  the  zero-order  Hankel  transrorm  of  the  Initial 
impulse  di-striiiution  I(r),   It  is  given  by 

00 

(2.6)  Kg")   =  y  I(r)J^(^)r  dr  . 

o 

Thus  to  compute  A,  T  and  \   for  a  given  value  of  r//gH  t 
we  merely  find  cf  f rom  (2,2)  and  Insert  it  into  (2,3)  -  {2»S).      Com- 
putationally it  is  much  easier  to  choose  a  value  of  <f   and  insert  it 
into  (2.2)  -  (2.5),  thus  determining  corresponding  values  of  r//gh  t, 
A,  T,  and  X    .   To  facilitate  this  computation  graphs  of  ^{d)    and  of 
/'•  \,  J  j  I ^)        are  given  in  Figures  1  and  21.   By  this  procedure,  we 
can  construct  a  graph  of  |a|  as  a  function  of  r//gh  t  (Figure  31). 
this  single  graph  gives  the  amplitude  of  the  waves,  except  for  the 
factor  I  yR^/o/g"  r  ,  as  a  function  of  r  for  any  fixed  value  of  t 
or  as  a  function  of  t  for  any  fixed  value  of  r. 

The  amplitude  A  tends  to  zero  as  tends  to  zero  if 

v^  t 
the  total  energy  of  the  wave  motion  is  finite.   This  energy  Q  is 

given  by 


00 

(2.7)  """P    i  ^ 


Q  =  ^  J   Ul(s)]''  tanh  sh  ds 


If  I(r)  is  discontinuous  then  Q  is  infinite.   This  is  the  case  if 
I(r)  is  a  non-zero  constant  within  a  circle  of  radius  R  and  zero  outside 
it.   V/e  will  consider  an  example  of  this  type  belbw.   In  the  other 
case  (of  initial  displacement)  the  total  energy  is  always  finite  pro- 
vided that  a  finite  volume  of  water  is  initially  displaced. 


-.•■;-i:)- 
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The  most  important  features  of  the  sunplitude  curve  seem 
to  be  the  location  and  magnitude  of  the  maxima  of  (aI.   These 
maxima  are  determined  by  (E«5)  and  therefore  depend  upon  the  ini- 
tial impulse  distribution.   However  a  simple  upper  bound  on  |a|  can 
be  found  for  any  initial  shape.   It  is 

(2.8)       |a|   <  1.28       if  R/h  <   2   . 

In  deriving  (2.8)  it  was  assumed  that  the  "effective  radius"  R  is 
defined  by 

CO 


(2.9)       I  ll^l  R^  =   r  r  ll(r)|  dr 


With  this  definition,  ttJI^Ir  is  just  the  total  initial  impulse. 

The  amplitude  A  is  zero  for  r  >  v/gh  t  and  is  continuous 
at  r  = v/gh  t  ,  though  its  derivative  has  a  jump  there.   Thus  the 
height  of  the  waves  increases  gradually  from  the  value  zero  in 
back  of  their  forward  boundary  r  =  ygh  t  . 

In  Table  I  the  maximum  amplitude  A    ,  the  group  velocity 

v„„„   of  the  maximum,  the  wave  length  X     and  period  T     at  the 
J^^x  '  °    max      ^       max 

maximum,  and  the  total  energy  Q  are  given  for  five  different  impulse 
distributions  I(r).   The  first  of  these  distributions  was  referred  to 
above  and  produces  a  motion  with  infinite  energy.   Consequently  the 
amplitude  has  infinitely  many  maxima  which  have  about  the  same  magni- 
tude and  occur  Increasingly  often  as  t  increases.   The  other  four 
distributions  all  lead  to  motions  with  finite  energy.   The  last  three 
have  no  secondary  maxima  while  the  second  has  infinitely  many,  the 
magnitudes  of  which  decrease  like  t    (for  fixed  r). 


1    ':  .> 
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TABLE 

I 

Kr) 
^o 

1,    r<R 
0,    r>R 

0,                  r>/2R 

r   2   -3/2 
[l+2(g)^] 

-y2rR-l 
e 

,-(§,2 

^MAX 

1»1 

0.79 

O.l+l 

0O5 

0,76 

V 

./gR 

0.33 

0.30 

0.35 

0.1+2 

0.38 

^MAX 
R 

2.8 

2.3 

3.1 

k'k 

3.6 

T 

•^MAX 

v/R/g 

k'l 

3.7 

k'3 

5.3 

I4..8 

Second- 
ary 
Maxima 

Infinitely 
majiy;    prac- 
tically 
constant 
amplitude . 

Infinitely 
many;    ampli- 
tudes  decrease 
like    t-2 

None 

None 

None 

Q/ttR/>"^i; 

CO 

> 

0.14-8 

0.18 

O.lli 

0.31 

i(s) 

^0 

|j,(Rs) 

^J2(/2Rs) 
s 

lj,2^-Rs//I 

r2 

^R'e'^^'^% 

e[l4R232]3/2 

-9- 

The  existence  of  secondary  maxima  and  their  rate  of  decrease  de- 
pend upon  the  smoothness  of  I(r),   If  I(r)  is  discontinuous  in- 
finitely many  maxima  of  practically  constant  amplitude  will  occur. 
If  I(r)  is  continuous  but  its  derivative  discontinuous  for  some 
r  =1=  0,  infinitely  many  maxima  will  occur  with  amplitudes  decreas- 
ing  like  t    (for  fixed  r).   If  I(r)  and  all  its  derivatives  are 
continuous  there  are  usually  only  a  finite  number  of  maxima,  and 
if  I(r)  is  a  decreasing  function  only  one  maximum  can  be  expected. 

A  graph  of  |A(r//gh  t)|  is  shown  in  Figure  31  for  the 
parabolic  impulse  distribution 

[1  -  i(§)2]  ,     r  <  /2  R  , 
(2.10)     I(r)   =  {      '^  ^ 

0  ,     r  >  Y^  R 

2 
described  in  the  second  coltomn  of  Table  I,  taking  R  =  —  h  ,   In 

Figures  l^Ia  and  L|.Ib  the  wave  height  >^(r,t)  is  shown  as  a  function 
of  t  for  two  different  fixed  values  of  r,  while  in  Figure  5l  ^  is 
shoivn  as  a  function  of  r  for  a  fixed  t. 

Similar  curves  for  other  values  of  r  or  t  can  be  ob- 
tained by  using  Figure  31  and  equation  (2,1).   To  obtain  correspond- 
ing curves  for  any  other  value  of  R/h  the  graph  of  A  for  this  vfclue 
of  R/h  must  be  constructed  from  (2.5).   To  facilitate  this  construc- 


/M)tahh  (f 


tion  the   graph  of  /      l^i  ( rf\ i^  Figure    21   should   be   used.      In 

')tanh 


calculating  A,  one  must  multiply  / ^{d)tar)h.  <S     by 


o  o 
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For  the  other  impulse  functions,  of  course,  one  may  employ 
the  same  procedure,  using  the  Hankel  transforms  1   given  on  the  bot- 
tom line  of  Table  I. 

3.   Initial  Elevation  or  Depression  of  the  Surface. 

The  height  of  the  surface  produced  by  an  initial  elevation 
or  depression  is  given  by 

(3.1)   Vj(r,t)  ^  Z£-.  B  cos  2M^  -   |)      for  r  »  R  . 

The  amplitude  B,  the  period  T  and  the  wavelength  X  are  functions  of 
,  while  h^  is  the  initial  elevation  at  the  origin.   R,  the 

yiH  t  '         f° 

effective  radius  of  the  initial  displacement,  is  defined  by  (3»5)« 
Both  X  and  T  are  given  by  the  same  equations  as  in  the  previous  case 
[equations  {2,3>)    and  (2.[^.)  respectively],   B,  which  is  zero  for 
r  >  ygK  t  ,  is  given  by 


<3-2)  ^  =  ^  2<H>/?^  for  r< /girt     . 

The  function  E(r-)  is  the  zero-order  Hankel  transform  of  the  ini- 
tial elevation  E(r)  and  is  given  by 

00 
(3.3)      S(|)  =  J      E(r)  J^i%)T   dr  . 


{:'■■:■  :  -v 


Ti.   .' 


■:.i-'- 


;3ut-- 
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Again  the  precise  form  of  the  function  B  depends  on  the 
exact  shape  of  the  initial  elevation,  but  again  a  simple  upper 
bound  exists  independently  of  this  shape.   We  find 

(3.I1.)     1b|  <  1.16   provided  R/h  <  1.6i|  , 

with  a  definition  of  R  analogous  to  (2.9 )  ,nairxely 

00 


(3.5)    I  |V|q|R^   =  t    |E(r)|r  dr 


With  this  definition,  Tr|Vi  JR  becomes  just  the  vol\arae  of  water  dis- 
placed by  the  given  initial  elevation  function  E(r). 

The  cjualitative  features  of  the  wave  motion  are  generally 
similar  to  those  described  in  the  previous  case.   A  significant 
difference  is  that  here  the  wave  amplitude  B  is  not  continuous  at 
r  =  y/gh  t  ,  so  that  the  waves  are  preceded  by  a  bore  or  discontinuous 
change  of  height  of  "dimensionless  amplitude" 

and  duration 


(3.6a)  T^^^^      =     l.lih^/V/Vyi       , 


where  R'  is  defined  by: 
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00 


(3.7)      I  \%\    r'^  =   1e(0)1   =   iJrE(r)  dr  1  . 

o 

The  sign  in  (3»6)  is  the  same  as  that  of  E(0).   Behind  the  bore 
the  amplitude,  after  perhaps  a  further,  continuous  increase, 
diminishes  gradually,  becoming  zero  at  r  =  co  .   Actually,  of 
course,  the  amplitude  is  continuous  at  r  =  /gh  t  and  it  is  not 
exactly  zero  for  r  >  /gh  t  .   The  bore  really  represents  a  very 
rapid  change  in  A  from  zero  to  the  value  given  by  (2.9) »   As  r 
increases,  the  bore  becomes  steeper,  approaching  a  discontinuity. 
The  bore  is  absent  if  the  average  initial  wave  height  is  zero* 
This  is  the  case  if  the  volume  of  fluid  above  the  undisturbed  water 
surface  is  equal  to  the  volume  of  fluid  which  has  been  removed  from 
beneath  the  undisturbed  surface,  i.e.,  if  E(0)  =  0  .   Since  this 
will  generally  be  the  case  in  explosion  problems,  no  bore  will 
usually  occur  in  these  problems. 

As  an  example  we  consider  the  initial  elevation  (or  de- 
pression) 

r  <  R   ■ 


A  graph  of  lB(r//gh  t)|  is  shown  in  Figure  3E  for  R  =  |-  h  .   In 
Figure  l|JE  the  i^rave  height  Vi(r,t)  is  shown  as  a  fvmction  of  t   for 
r  =  20  km  .   In  Figure  5E,  h  is  shown  as  a  function  of  r  for 
t  =  7.5  min.   Similar  curves  for  other  values  of  r,  t,  or  R/h,  or 
for  other  initial  shapes,  can  be  constructed  from  Figure  2E  and  3E 
and  equations  (3*1)  and  (3*2)  just  as  in  the  impulse  case. 
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For  the  example  of  (3«8)  we  find  the  following  results; 


TABLE  II 


Maximum  Amplitude  B„ 

max 

Velocity  of  Maximum  V 

''  max 

Wavelength  at  Maximum  X 

Period  at  Maximum  T 

max 

Amplitude  of  bore  B, 
Velocity  of  bore   V 


max 


Duration  of  bore   T 


bore 


»82 

.37/iR 
3.5R 
Il..6/R7g 

-  2h 


l.ilh^/^r^/V/g 


The  maximum  amplitude  listed  above  occurs  at  the  first 
maximum  of  1b|.   There  are  infinitely  many  subsequent  maxima  the 
amplitudes  of  which  decrease  like  r//^  t  .   The  amplitude  and 
velocity  of  the  second  maximum  are 


(3.9) 


Second  Maximum   B    ^  =  .[1.9 

max  c 

Velocity  of  Second  Maximum  V   p  =  .2^/gR 


Returning  again  to  the  general  case,  we  consider  the 
limiting  form  of  our  result  when  h  becomes  infinite.   From  (2.2) 
we  then  find  that  cf/h  =  gt  /l\.v      .   Inserting  this  into  (3«1),  (2.3) 
-  (2.1}.)  and  (3,2)  yields 


(3.10) 


i^(r,t) 


^o^ 


f  B  cos  0- 


-Ih- 


(3.12)  X      =  ^ 

(3.13)  T     =    ^ 


The  amplitude  of  the  bore  given  by  (3 •6)  tends  to  zero.   The  re- 
sults (3,10)  -  (3.13)  include  the  special  cases  of  parabololdal  and 
ellipsoidal  depressions  treated  by  Poisson  [3]« 

l\..        Waves  Produced  by  Explosions* 

In  order  to  apply  the  preceding  results  to  the  production 
of  waves  by  an  explosion,  we  must  determine  the  initial  Impulse  dis- 
tribution in  the  first  case,  and  the  Initial  displacement  of  the 
surface  in  the  second  case.   The  initial  impulse  distribution,  for 
an  above  surface  burst,  can  be  computed  from  the  pressure  pulse 
(shock  wave)  emitted  by  the  burst.   However,  the  initial  displace- 
ment cannot  be  so  readily  deduced.   Nevertheless  estimates  of  this 
displacement  can  be  made.   Let  us  begin  with  the  Impulse  case. 
A.   Initial  Impulse 

Suppose  that  an  explosion  occurs  at  a  height  H  above  the 
water  surface.   Then  the  distance  from  the  explosion  to  a  point  on 


the  surface  Is  y/Ti  +r   .   Here  r  denotes  radial  distance  from  a  point 
on  the  surface  just  below  the  explosion  ("surface  zero").   Now  let 

jl  WE   +r  Jbe  the  impulse  per  unit  area  of  the  incident  shock  at 

7"2  ? 
distance  /H  +r  from  the  explosion.   Then,  neglecting  the  effect  of 

shock  reflection  at  the  water  surface,  we  have 


■1\- 


;>.;.:      .V.:.      i-i\. 
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^      ^  Ji(H) 

Actually  I(r)  will  be  greater  than  (I|..l)  indicates  due  to  reflection. 
This  effect  will  be  greatest  for  small  values  of  r  and  decreases 
rapidly  as  r  increases. 

Let  us  assume  that  ^  l/h  +r  I   is  inversely  proportional 
to  (/h  +r  p.  (The  peak  pressure  in  the  Taylor  point  blast  varies  in 
this  vjay,  although  the  impulse  does  not.)   Then  (i+.l)  becomes 


(1^*2)        I(r)  =J{E)    [1   +  (§)2]-^/2 


This  is  just  the  case  treated  in  column  four  of  Table  I  provided  we 
set  R  =  s/z   H  *   Thus  in  this  case  the  wave  pattern  contains  a  single 
maximum  which  moves  outward  with  the  velocity  V    =  05  v^  /s^  — 
»l\.^/W-    i   The  wavelength  at  this  maximum  is  X    =  3*\/2   H  =  l+.l^H 
and  the  period  T^^^^  =  i|.*3  ii/2  v/H/g  =  5*\/Wg    •   The  wave  height  at  the 
maximum  is 


The  energy  in  the  wave  motion  is 

ik^k)  Q  =  1. Sir  V2x/H<9  2(h)/ -1  =  .67yHi^2(H)/3-^   . 

Equation  {k»5)    shows  that  the  maximum  wave  height  depends  upon  yH^(H) 


/ 
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and  therefore  it  increases  as  H  decreases*   The  total  energy  Q  in 
the  wave  motion  also  increases  as  H  decreases.   Of  course  the  above 
results  are  invalid  at  H  =  0  and  therefore  also  for  small  values  of 
H  since  they  become  singular.   This  is  due  to  the  improper  choice 

Of  J'fyS^)  . 

More  reliable  results  can  be  obtained  by  using  a  more 
accurate  choice  of  I(r). 

B*   Initial  Displacement 

The  potential  energy  \-J   contained  in  the  water  initially 
is 

9° 

(i]..5)        W  =  2iT^g  J   E^(r)rdr  . 

o    -^ 

Let  us   first    assxame  that  E  =  Vi     for  r  <  R  and  E  =  0  for  r  >  R  as 
in   (3.8)-      Then  W  =  rr  p  gh^R     and    (3.1)   becomes 

According  to  Table  II  Ib  _|  =  .82  so  we  have 
°  max 

1/2 
(U.7)         |V^(r,t)|  <  .82(1^)     .  i  . 

Other  properties  of  the  wave  motion  are  given  in  Table  II  and  other 
parts  of  section  3«  Of  course  V/  must  be  known  or  estimated  in  some 
way  in  order  to  use  these  formulae  to  predict  vjave  heights.  On  the 
other  hand,  from  observations  of  wave  height  W  can  be  computed. 
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Now  suppose  that  the  Initial  displacement  is  arbitrary, 

and  not  necessarily  that  just  considered.   Let  V  denote  the  volume 

between  the  initial  surface  and  the  undisturbed  surface.   Then  from 

,2 


(3*5)  we  have  njh  [r  =  V  ,   Thus  (3.1)  bee 


omes 


(i^„8)  yr,t)  -  ^  .  I  cos  271  (|  -  I) 


According  to  i3»k)    |b|  <  1.16  provided  that  R/h  <  1.61;.   Thus  {3»k.) 
yields 


(1^.9)        lV^'*)l  <  ^^W^  •  F        (R/h  <  1.61^) 


Again  (i+.9)  can  be  used  to  estimate  wave  height  only  if  V  and  R  are 
known  or  estimated. 

In  case  some  water  is  removed  from  the  ocean  by  the  ex- 
plosion, a  negative  bore  will  occur,  preceding  the  other  waves.   This 
may  happen  if  water  is  evaporated  or  thrown  into  the  air  and  retained 
there  in  the  form  of  droplets.   The  bore  travels  with  the  velocity 
y^  .   The  duration  of  the  negative  phase  behind  the  bore,  T^qj,^  is 
given  by  (3.6a).   For  the  explosion  of  Figure  3E,  at  r  =  20  Km, 
T,     =  5l  sees.   This  can  also  be  seen  from  Figure  l+E.   At  r  =  2500Km, 

Tv    is  five  times  as  great.   The  maximum  follows  the  bore  by  the 
bore  ° 

delay  time  — -—  ^ 


^S^  -  ^max 


'•  i 


7  v.:::-  -f 
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Appendlx.   Mathematical  Derivation  of  the  Results. 

Consider  an  incompressible  fluid  of  density  f    and 
depth  h  over  a  rigid  plane  bottom.   Suppose  that  at  t  =  0  a  cylin- 
drically  sjnTimetric  distribution  of  impulse  I(r)  is  applied  to  the 
surface,  which  is  initially  horizontal  and  at  rest*   Then  the  fluid 
will  begin  to  move  and  its  surface  will  no  longer  be  horizontal. 
Let  Vi(r,t)  denote  the  height  of  the  surface  above  its  initial  posi- 
tion.  Then  according  to  the  linear  theory  of  surface  waves, 
h(r,t)  is  given  by 

00 

(A.l)    h(r,t)  =  -  -i-  lim     /   s^/^I(s )ytanh  sh 
^         PJH    y->0"  J 


/Vi  y 


sech  sh  cosh  s(y+h)sin[/gs  tanh  sh  t]  J  (rs)ds  « 


This  integral  representation  of  the  wave  height  Vi  was 
obtained  by  applying  the  Hankel  transform  to  the  Laplace  equation 
for  the  potential  function  which  describes  the  flow  (see,  e.g. 
Stoker  [15]  or  Sneddon  [16]).   Of  course  the  appropriate  initial 
and  boundary  conditions  were  also  utilized.   Then  from  the  potential 
function  the  wave  height  was  calculated. 

The  integral  in  the  above  expression  for  h  can  be  evalu- 
ated asymptotically  for  large  r  and  t.  To  do  this  we  first  replace 
the  Bessel  function  J  (rs)  by  its  asymptotic  expression  in  terms  of 
a  cosine,  which  is  valid  when  the  argument  is  large.   Then  we  assume 
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that  the  ratio  t/r  is  fixed"  ,  so  that  the  resulting  integral  con- 
tains just  one  large  paraineter,  either  t  or  r.   Finally  we  apply 
Kelvin's  stationary  phase  formula  [11]  assximing  that  t  or  r  is  large. 
In  this  vray  v;e  obtain  equation  (2.1)  together  with  the  related  equa- 
tions (2.2)  -  (2.5). 

The  bound  on  |a|  given  by  equation  (2.7)  follows  from  (2.5) 
and  certain  properties  of  the  Hankel  transform.   To  deduce  it  we  first 
note  that  the  function  fjicS)    =  J^^^^^^f^    has,  for  large  cf,    the 
asymptotic  expansion 

(A. 2)  \l/j(crO  ~  y2^  (1  -  8  cf^e"^°'+  ...)  . 

This  indicates  that  ^j{(^)   may  be  less  than  its  asymptotic  value  s/25'  . 
Now  \|/-[-(d)  is  bounded  at  di   =  0^  is  continuous  for  all  non-negative  o*  , 
and  has  the  pov;er  series  expansion 


(A. 3)  \l»j(c<)   =  1  +  i  c^  +  ... 


From  this  expansion  it  appears  that  ^j{c^)    <   /2tf  for  c<  >  0.53  (cf. 
Figure  21).   Using  this  inequality  in  (2.5)  now  yields 


(A.li.)      Ul  <  — jZI^^^  cf3/2  |j(rf),       o'>0.53  . 

o 


In  the  case  of  water  of  infinite  depth  it  is  usual  to  assume 
that  t  /r  is  constant  rather  than  t/r,  as  in  the  present  case 
(cf.  Stoker  [15]  or  Sneddon  [16]). 
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If  we    set    s   =  rf/h  this  becomes 

(A.5)  lA|     <        .   ^^    .   ;p       B^^^\i{s)\ 

Nov;  -  sl(s)  is  the  first  order  Hankel  transform  of  the  derivative 
I'(r),   This  follows  from  an  application  of  the  Hankel  inversion 
formiila  and  differentiation  under  the  integral  sign.   Therefore  (A«5) 
may  be  written  as 

(A. 6)         |a1  <  ^  ^{^.y„  s^/2{  (  rl«(r)j.(rs)drl. 

[iJr1/2       J  1 

Suppose  |l(r)j  decreases  monotonically  to  zero,  as  must  be  the  case 
in  an  explosion  problem.   Then  v;e  may  evaluate  the  Integral  in  (A. 6) 
by  employing  a  simple  estimate  for  the  Bessel  function  J,  ,  which 
yields 

CO  00 

(A. 7)     |a!  <  --2 [  r^/2j,(^^^^  ^  _1__   Tpl/Sj^j,)^^  ^ 

If  we  now  set  K  =  v^R  we  may  write 

K  00 

(A. 8)    |a|  <  — i—    f  r"^/2j(^)^^  +  _1_    fr*^/2j(^j^^ 

~  lyfrR  X  lytrR  ,i 

o  o  o  K 

K  CD 

<  -i-.    fr'^/^dr  +  \/^  rrl(r)dr     =     1.28      . 

This   is  formula   (2.7)    for  6-  >  0.53 .      If  RA  is   so   small   that   the   first 


, ■»•  i  .i''.» 
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extremvun  of  <3I(r')  occurs  at  a  value  of  tf'  greater  than  0.53,  (A. 8) 
holds  for  all  d.      Thus  the  bound  in  (2.7)  has  been  deduced. 

Let  us  now  consider  the  case  of  aja  initial  elevation  or 
depression  of  the  surface.   In  this  case  the  wave  height  is  given  by 

oo 
(A. 9)   h(r,t)  =  „  vo  ~  T     sE(s)  sech  sh  cosh  3{y+h)cos[/gs  tanh  sh  t] 

°  jQ(r3)ds. 


The  asymptotic  formulas  (3.1),  (3.2)  are  derived  from  (A. 9)  by  station- 
ary phase  in  exactly  the  same  way  as  for  the  Impulse  case. 

Owing  to  the  possible  existence  of  a  bore,  however,  the 
upper  bound  on  |b|  must  be  derived  somewhat  differently.  Since  the 
factor  E(t-)  in  (3*2)  does  not  necessarily  vanish  at  o"  =  0  [as  did  the 
corresponding  factor  cfl(r-)  in  (2.5)1,  an  upper  bound  on  the  comple- 
mentary factor  \j/p(o')  ~  y^fir^fT  must  be  found  which  is  valid  for  all 
positive  cf  ,  not  merely  for  cf  sufficiently  large.  Again  this  bound 
is  suggested  by  the  series  expansions  of  v|/„(o')  about  the  origin: 


(A. 10)  *£(c<)   =  1  +  ^  c^  +  ... 


and  near  infinity: 


(A. 11)  ^^{d)      =    ^2   cf[l   -   8cr^e"^^  +    ...] 


It   is 


(A. 12)      \)/g((J)   < 
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J2  for  0  <  d  <  1   , 

6/1  for  d  >  1   , 


-  cf.  Figure  2E.   Thus  from  formula  (3.2) 

(A. 13)  |B|      <      jl^      I8(e>1  ior4<X      , 

tl^E     <'IE(|)|        fox.c^>l        .      . 

If  monotonlcity  of  the  Initial  elevation  E(r)  is  assumed,  the  first 
line  yields  at  once 

(A.ll^)      1b1   <   [^[p^  |e(0)  1  =  0.71  R/h     for  d  <  1   , 

while  the  second  can  be  broken  up  as  in  (A, 6)  -  (A, 8): 

oo 
(A. 15)         1b1      <     -|^i^s|E(s)|g     =     1#^    I    I    rE'(r)J^(rs)dr| 

00  CB 

-^^     j       .K.(r)dr=M|    j     E(r)dr 
K     °  oo^       ^ 

t°      g  1o      K 


00 

l^^  +  ^^ll     f      rE(r)dr=l,l6  for   c^  >  1      . 

^        o 

Since  the  bound  of  (A,l$)  is  larger  than  that  of  (A.lii.)  whenever  R/h 

<  I.6I4.,  formula  (3«i|)  is  established. 
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Pormula  (2.7)  for  the  energy  of  the  water  motion  in  the 
Impulse  case  should  be  just  as  well-known  and  standard  as  its  eleva- 
tion counterpart  (JH..5),  but  this  seems  somehow  not  to  be  the  case 
(cf.  Unoki  and  Nakano  [1[|.]).  Its  derivation  affords  an  instructive 
example  of  the  direct  use  of  the  exact  integral  representation  of 
the  motion.  At  t  =  0  all  of  the  energy  of  the  motion  is  kinetic. 
Hence  the  total  energy  Q  is  given  by  a  Dirichlet  integral 

o  00 

(A. 16)     Q  =   1  J  y^  2Try^'"[ir  "^  ly^^  <^^  ^^ 

-h   o 

of  the  velocity  potential  $(r,  y,  0)  over  the  entire  volume  occupied 
by  the  water.   By  Green's  formula  and  the  bottom  condition,  this 
simplifies  to 

CO 

(A. 17)      ^     "^     ^f    f   rf(r,0,0)|  (r,0,0)dr. 

o 

But  in  the  course  of  deriving  (A.l)  one  finds  that 

f 

(A. 18)     $(r,y,o)  =  -  i  J  sl(s)  sech  sh  cosh  s(y+h) jQ(rs)ds. 


r 


Hence 


(A. 19)  |(r,0,0)   =   -  i  I(r) 


and 


-2k' 

CX) 

(A. 20)      I  (r,0,0)  =  -  •!  J  s^i(s)tanh  sh  J^(rs)ds. 

Substituting  into  (A.l?)  and  interchanging  the  order  of  integration, 
one  finds  that 

oo  00 

(A. 21)   Q  =  ?  f  s^I(s)tanh  sh  f  rl(r ) jQ(rs)dr  ds 


o 

00 


=  5  /  [sl(s)]^tanh  sh  d 


s  , 


as  in  formula  (2.7)   . 
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